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ABSTRACT 

From a sample of ~ 50, 000 early- type galaxies from the Sloan Digital Sky Survey 
(SDSS), we measured the traditional Fundamental Plane in the g, r, i and z bands. 
We then replaced luminosity with stellar mass, and measured the "stellar mass" Fun- 
damental Plane. The Fundamental Plane, R oc a a / I B , steepens slightly as one moves 
from shorter to longer wavelengths: the orthogonal fit has slope a = 1.40 in the g band 
and 1.47 in z, with a statistical random error of ~ 0.02. However, systematic effects 
can produce larger uncertainties, of order ~ 0.05. The Fundamental Plane is thinner 
at longer wavelengths: it has an intrinsic scatter of 0.062 dex in g and 0.054 dex in 
z. We have clear evidence that the scatter is larger at small galaxy sizes/masses; at 
large masses measurement errors account for essentially all of the observed scatter 
(about 0.04 dex), suggesting that the Plane is rather thin for the very massive galax- 
ies. The Fundamental Plane steepens further when luminosity is replaced with stellar 
mass, to 1.54 or 1.63 when stellar masses are estimated from broad-band colors or 
from spectra, respectively. The intrinsic scatter also reduces further, to 0.048 dex on 
average. Since color and stellar mass-to-light ratio are closely related, this explains 
why color can be thought of as the fourth Fundamental Plane parameter. However, 
the slope of the stellar mass Fundamental Plane remains shallower than the value of 
2 associated with the virial theorem. This is because the ratio of dynamical to stellar 
mass increases at large masses: M^/M* oc M^ 7±0 ' 01 . This scaling is the edge-on 
projection of the stellar mass ft-space. The face-on view suggests that there is an up- 
per limit to the stellar density for a given dynamical mass, and this decreases at large 
masses: M*/Rl oc M d yl 3 . All these trends can be used to constrain early-type galaxy 
formation models. 

We also study how the estimated coefficients a and B of the Plane are affected 
by other selection effects, whether in apparent or absolute quantities. For example, 
if low luminosity objects are missing from the sample, and one does not account for 
this, then a and B are both biased low from their true values. If objects with small 
velocity dispersions are missing, then a is biased high, although this matters more for 
the orthogonal than the direct-fitted quantities. These biases are seen in Fundamental 
Planes which have no intrinsic curvature, so the observation that a and B scale with 
L and a is not, by itself, evidence that the Plane is warped. On the other hand, we 
show that the Plane appears to curve sharply downwards at the small-size/mass end, 
and more gradually downwards as one moves towards larger sizes/masses. Whereas 
the drop at small sizes is real, most of the latter effect is due to correlated errors. 
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1 INTRODUCTION 
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The Fundamental Plane has been a useful diagnostic of 
galaxy distances and galaxy evolution (e.g. Dressier et al. 
1987; Djorgovski & Davis 1987; J0rgensen et al. 1996; Pahre 
et al. 1998; Bernardi et al. 2003; J0rgensen et al. 2007). If 
galaxies are virialized, then one expects observable signa- 
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tures of the balance between potential and kinetic energies. 
For example, 



2 GMdyn Mdyn L Mdyn T 



(1) 



where a is a velocity dispersion, / is a surface brightness, 
R is a scale, and Mdyn/L is the mass to light ratio. This 
suggests that the observed line-of-sight velocity dispersion 
a\ and surface brightness I e = L/(2ttRI) at or within some 
fiducial radius R e should be correlated with one another. (In 
what follows, we will use the subscript e to denote quanti- 
ties estimated from deVaucouleur's (1948) fits to the surface 
brightness profiles.) 

Velocity dispersion and surface brightness are distance- 
independent observables, whereas the size is not, so it is 
common to write 



Re 



(2) 



and to find that pair (a, B) for which the scatter either in 
the R e direction, or in the direction orthogonal to the Plane, 
is minimized. If the resulting values of (a, B) differ from the 
virial scalings (2, 1), then the Fundamental Plane is said to 
be 'tilted'. The tilt of the Fundamental Plane is interpreted 
as evidence that the mass-to-light ratio (Mdyn/L) depends 
on some combination of the observables (R e , I e ,cr e ) - the 
prejudice that galaxies are virialized sets what this com- 
bination must be: (Mdyn/L) oc G 1 e ~ a i^f -1 ". In low redshift 
(z rsj 0.1) samples, a ^ 1.5 and B ~ 0.8, so the implied 
Mdyn/L is expected to increase slightly with Md yn or L. 
J0rgensen et al. (2007) find that a = 0.6, B = 0.7 and 
Mdyn/L oc M - 54 at z ~ 0.9. 

Contributions to the tilt of the Fundamental Plane can 
be further examined by expressing the mass-to-light ratio in 
terms of the dynamical mass, total mass, stellar mass, and 
broadband luminosity: 



M, 



dyn 



M, 



dyn 



Mtot 
M» 



T~ ) (3) 



where Mtot is the sum of dark matter and baryonic mass, 
M* is the stellar mass, and L is the observed broadband 
luminosity. If any of these terms varies as a function of Mdyn 
or L, the Plane will be tilted. 

The assumption that Md y n/M to t is constant is the as- 
sumption of homology. The validity of this assumption is 
supported by detailed kinematic modeling of galaxies ob- 
served with SAURON integral-field spectroscopy. The com- 
bination of 2-D photometric and spectroscopic data allowed 
Cappellari et al. (2007) to account for differences in e.g. 
p(r), cr(r), and V(r). They found that dynamical mass is 
a robust tracer of total mass. More recent work by Bolton 
et al. (2008) suggests that the deVaucouleur-based quantity 
ReCTe is linearly proportional to the mass estimated from the 
strong gravitational lensing effect, further reinforcing the ho- 
mology assumption that Md yn oc Mtot- 

For a given IMF, the stellar mass to light ratio, (M*/L), 
depends on the age and metallicity of the stellar population 
as well as on wavelength (e. g., Tinsley 1978; Worthey 1994). 
There is now growing evidence for correlations between mass 
and metallicity (e.g., Trager et al. 2000; Nelan et al. 2005; 
Thomas et al. 2005), and between age and velocity disper- 
sion (e.g., Bernardi et al. 2005). Hence, stellar population 
models predict that (M*/L) will depend on total mass and 



luminosity, plausibly producing a tilt that will depend on 
waveband. While there is evidence that the tilt does indeed 
depend on waveband, the dependence is weak (Pahre et al. 
1998; Bernardi et al. 2003; La Barbera et al. 2008), and so 
this effect alone cannot explain all of the tilt. 

There is little discussion in the literature of the pos- 
sibility that the IMF changes along the sequence in such a 
way as to make the dependence on M*/L weak. If there is no 
such effect, then stellar population effects and non- homology 
are weak, so the dominating contribution to the tilt of the 
Fundamental Plane is due to the systematic variation of 
Mtot I 'M* with mass. In a hierarchical scenario of galaxy for- 
mation, the relative distribution of stars and dark matter in 
a galaxy depend on the roles of dissipational and dissipation- 
less merging (e.g. Bender et al. 1992). Hydrodynamical sim- 
ulations of galaxy mergers (Robertson et al. 2006; Hopkins 
et al. 2008) suggest that the fractional gas content of merg- 
ing galaxies determines the fundamental Plane of their rem- 
nants. Further dissipationless mergers preserve the Funda- 
mental Plane, but not its projections (e.g., Boylan-Kolchin 
et al. 2005). Observations of brightest cluster galaxies by 
Bernardi et al. (2007) support these conclusions. 

The dissipational content in mergers sets the effective 
mass-to-light ratio of the merger remnant (Robertson et al. 
2006; Hopkins et al. 2008). If the dissipational content varies 
with mass (in spirals the gas fraction decreases with mass, 
e.g. Bell & deJong 2000), this would cause a systematic de- 
pendence of the mass-to-light ratio with mass. If this ratio 
varies as a power law with mass (M/L oc M a ), then the Fun- 
damental Plane would be tilted relative to the virial scaling. 
If the mass-to-light ratio varies in a more complicated man- 
ner, the Fundamental Plane could be warped, or have a more 
complicated shape. 

Therefore, in the present work, we separate out the con- 
tribution from these effects by re- writing equation (1) as: 



dyn 



L 



R. 



(4) 



We do not consider variation of (Mdyn/M to t), since two 
rather different methods (Cappellari et al. 2007; Bolton et al. 
2008) suggest that this ratio is a constant across the early- 
type population. (This is also why we do not include effects 
associated with fitting Sersic's 1968 generalization of the 
deVaucouleur profile to the images, even though there is a 
well-developed literature on this subject.) If the stellar popu- 
lation models which one uses to estimate M*/L are accurate, 
then multiplying the surface brightness by M*/L should 
eliminate most of the dependence on waveband. (There is a 
small remaining dependence which arises from the fact that 
the half-light radii are slightly but systematically smaller 
in redder bands, e.g. Hyde & Bernardi 2009). The remain- 
ing term (Md yn /M*) is the ratio of the dynamical to stellar 
mass, which we will use as our proxy for structural differ- 
ences. This is interesting because Md yn /M* is not expected 
to be very much greater than unity - recent work suggests 
that Mdyn / M* ~ 1/0.7 (Capellari et al. 2007) - how it scales 
with mass is a quantity of great interest to galaxy formation 
models (e.g. Bower et al. 2006; DeLucia et al. 2006; Hopkins 
et al. 2008), and previous work in the SDSS suggests that 
Mdyn / M* is not constant across the early- type population 
(Padmanabhan et al. 2004; Gallazzi et al. 2006) 

Section 2 discusses how we select our sample and shows 



the result of fitting for the coefficients (a, B) of the tradi- 
tional Fundamental Plane in the SDSS g, r, i, and z bands. 
(Appendix A shows that our sample probably contains some 
non- early- type galaxies but they are too few to affect our 
conclusions.) We also study the possibility of detecting if 
the Plane is warped, or merely thick. This is prompted in 
part by recent work showing that the parameters of the fit- 
ted Fundamental Plane appear to depend on the luminosity 
and velocity dispersion range of the sample (D'Onofrio et 
al. 2008; Nigoche-Netro et al. 2009). We show that such de- 
pendences exist in thick Planes that are not warped - they 
are a consequence of not accounting for selection effects. We 
then discuss more reliable estimates of curvature along the 
Plane. 

Section 3 describes our estimates of M*/L and the result 
of fitting for the coefficients (a, (3) in 

R e oc a" YTj (5) 

where E = (M*/L)I. The difference between (a, (3) and 
(a, B) is due to stellar population effects; the difference be- 
tween (a,/3) and the virial scalings (2, 1) encodes informa- 
tion about M dyn /M*. If (a, /3) < (2, 1) then (M dyn /M*) in- 
creases with Mdyn or M* for the same reason that (a, B) < 
(2, 1) implies (Mdyn/ L) increases with Mdyn or L. 

A complementary analysis of the Fundamental Plane 
variables was introduced by Bender, Burstein & Faber 
(1992). This construction, which they termed ft-space, was 
criticized as being an "obfuscation, not a simplification" , by 
Pahre et al. (1998), primarily on the grounds that R e and 
I e depend on waveband. However, by replacing luminosity 
with stellar mass, what we will call ft* -space, most of the 
wavelength dependence is removed; only the weak depen- 
dence of R e on waveband remains. So it is interesting to 
examine how early- types are distributed in ft* -space. Sec- 
tion 4 presents the first analysis of ft* -space. The scaling 
of M* I Mdyn with Mdyn referred above is the edge-on view 
of the ft* -Plane. The face-on view provides a relation be- 
tween the maximum stellar density and dynamical mass of 
early- type galaxies. 

A final section summarizes our findings and discusses 
some implications. 



2 FP: THE TRADITIONAL FUNDAMENTAL 
PLANE 

We use the sample of about 50,000 early- type galaxies as- 
sembled by Hyde & Bernardi (2008). The sample is based 
on the SDSS-DR4, but with photometric and spectroscopic 
parameters updated from the SDSS-DR6 database. Briefly, 
objects have deVaucouleur magnitudes 14.5 < m r < 17.5, 
f racDev = 1 in both the g- and r-bands, and axis ratios 
b/a > 0.6. Appendix A shows that these cuts almost cer- 
tainly do not yield a pure early- type galaxy sample (e.g., 
the face-on analogues of the objects with fracDev=l and 
b/a < 0.6 are still in our sample). However, the non-early 
types in our sample are too few to affect our conclusions. 

Photometric parameters for the best-fit deVaucouleur 
surface brightness profiles are "corrected" for known prob- 
lems which arise from the SDSS sky-subtraction algorithm. 
In addition, there are some systematic differences between 
the velocity dispersions output by SDSS-DR6 and the 
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IDLspec2d reduction (Hyde & Bernardi 2009). Our veloc- 
ity dispersion estimates are simply the average of the two 
reductions (Section 2.1.2 discusses the dependence of the 
Fundamental Plane parameters on systematics) . We also 
compute the velocity dispersion for objects for which the 
SDSS pipeline does not estimate a (due to low S/N or the 
presence of weak emission lines, i.e. the status flag not- 
equal to 4). We select galaxies with velocity dispersions 
60 < a < 400 km s _1 . The velocity dispersions are then 
corrected to R e /8. Stellar mass estimates, from Gallazzi et 
al. (2005), are available for all these objects. Since these are 
actually stellar mass to light ratios multiplied by a lumi- 
nosity, and we have corrected the luminosities for the sky- 
subtraction problems whereas Gallazzi et al. did not, we 
apply a correction to their stellar masses to account for this 
effect. Where necessary, distances were computed from red- 
shifts assuming a Hubble constant of 70 km s— 1 Mpc -1 in 
a flat ACDM model with Q = 0.3. 

2.1 Fitting the Plane 

We fit the Plane as follows. We begin by writing the Funda- 
mental Plane as 

^4^) =aloe "(^ +b ^ + c > (6) 

where c = (log 10 R) — a (log 10 a) — b (/x e ), and \i e is the mean 
surface brightness within the half light radius, defined ex- 
plicitly as follows: 



fie = -2.5 log 10 (J e ) = -2.5 log 10 




= m + 51og 10 (r e ) + 2.51og 10 (27r) - 101og 10 (l + z)(7) 

where m is the evolution, reddening, and k-corrected appar- 
ent magnitude (refer to Hyde & Bernardi 2009 for details 
on magnitudes), and r e is the angular size in arcseconds. 
In this form, the virial scaling would follow (a, b) = (2, 0.4) 
because of the —2.5 term in the definition of fi e - Unless we 
say so explicitly, we work in logarithmic units. Therefore, 
following Bernardi et al. (2003), we will sometimes abuse 
notation by using R to denote log 10 (i? e /kpc), V to denote 
log 10 (cr/km s _1 ), and / to denote fi e . Thus, for example, 
eR,ey and ej denote the measurement errors on logi? e , 
log a, and ji e - 

The shape of the Fundamental Plane is determined by 
estimating a and b. This is done either by minimizing resid- 
uals in the R e direction, or in the direction orthogonal to 
the fit. In general the 'direct' and 'orthogonal' fit param- 
eters are different combinations of the mean values of and 
covariances between the variables log 10 R, log 10 a and fi e . 

In practice, naive estimation of these means and covari- 
ances (e.g. simply summing over the data without including 
other weight terms) may lead to biases induced by measure- 
ment errors (these usually affect the covariances) or by selec- 
tion effects (which bias the means and the covariances). The 
effects of both must be accounted-for to estimate the intrin- 
sic shape parameters a and b (e.g. Saglia et al. 2001). This is 
especially important when the FP is determined for galaxies 
in a magnitude limited sample (Bernardi et al. 2003). We 
do this following methods described in Sheth & Bernardi 
(2009). Analytic expressions which quantify the bias in the 
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Fundamental Plane due to the magnitude limit of the survey 
may be found there. 

Briefly, the effect of the magnitude limit is removed by 
weighting each object by V~^(L), the inverse of the vol- 
ume over which it could have been observed, and measure- 
ment errors are subtracted in quadrature from the measured 
covariances. The mean values of magnitude, size, surface- 
brightness, and velocity dispersion in the r-band are 

(M r /mags) = -20.99, (log R e /kpc) = 0.39, 

(fie) = 19.53, (log a/km s" 1 ) = 2.19. (8) 



With typical measurement errors 

/ <e?> 
£ (eie R ) 

\ (eie v ) 

( 0.0542 0.0162 
0.0162 0.0049 
V 



(eie R ) 

{el) 
{e R e v ) 




(9) 



we find that the intrinsic covariance matrix (i.e. corrected 
for measurement errors and the magnitude limited selection 
effect) is 

(Cu Cir Civ \ 
Cir Crr Crv I 
Civ Crv Cvv ) 

I 0.2947 0.0782 -0.0096 \ 

0.0782 0.0552 0.0189 . 
\ -0.0096 0.0204 0.0187 / 



(10) 



From this covariance matrix we obtain 

adir = 1.170 and 6dir = 0.303 
a or t = 1.434 and b or t = 0.315. 



(ii) 

(12) 



The observed rms scatter about the direct fit is 0.107, of 
which 0.096 is intrinsic. These quantities for the orthogo- 
nal fit are 0.066 and 0.058 respectively. The uncertainty on 
the coefficients a and b is dominated by systematic effects 
more than random errors, as described here below. Typical 
uncertainties on the coefficients due to random errors are 
5a < 0.02 and 5b < 0.01. Sytematics errors give 5a ~ 0.05 
and 5b ~ 0.02. 



2.1.1 Dependence of 5a and 5b on sample size 

We analyze how the uncertainty on a and b depend on sam- 
ple size using bootstrap resampling. Figures 1 and 2 show 
how 5a and 5b depend on sample size. They were obtained 
by dividing the total sample up into smaller subsamples, 
each containing iVg a i galaxies. For a given iV ga i, we do not 
include any galaxy more than once, so there are a total of 
N = N to t/N g& \ subsamples. We then fit the FP in each sub- 
sample, and computed the mean and rms values of a, b and 
intrinsic scatter. Figure 1 shows that the precision increases 
as the number of objects in the sample increases, as one 
might expect. However, Figure 2 shows that sample sizes 
smaller than about 300 tend to result in underestimates of 
the intrinsic scatter, in good agreement with La Barbera et 
al. (2000). For a sample of ~ 50,000 galaxies the statisti- 
cal random error 5a and 5b are quite small: 5a < 0.02 and 
5b < 0.01. 
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Figure 1. Precision of the FP coefficients a, b and scatter as the 
number of galaxies used to fit the FP changes. The quantities 
5a, 5b and 5FPrms show the standard deviation values of a, b 
and the scatter around the FP, obtained from fitting the FP to 
N = Nt t/N ga i subsamples of the whole sample (i.e. each galaxy 
was chosen only once). Triangles and circles show results for the 
direct and orthogonal fits, respectively. For clarity, the triangles 
have been slightly shifted to the right. 




Figure 2. As for Figure 1, but now for the dependence of the 
estimated intrinsic scatter around the FP on sample size. Error 
bars show the (5FPrms values from Figure 1. 
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Figure 3. The Fundamental Plane in the g, r, i and z bands. Typical uncertainties on the coefficients due to random errors are 5a ~ 0.02 
and 5b ~ 0.01. Sytematics errors give 5a ~ 0.05 and 5b ~ 0.02. 



Dependence of 5a and 5b on systematics 

When working with a large galaxy sample, systematic effects 
may be more important than random errors. The correc- 
tion applied to the magnitudes, sizes and stellar masses for 
known problems which arise from the SDSS sky-subtraction 
algorithm (Hyde Sz Bernardi 2009), results in small changes 
to a and b - the variation is less than 0.01 for both coeffi- 
cients. 

However, recall that our estimate of the velocity disper- 
sion is the average of the DR6 and IDLspec2d values. If we 
only use the DR6 values, we find (a d irect-DR6, frdirect-DR6) = 
(1.189,0.303) with intrinsic scatter 0.108 dex, and 
(a or th-DR6, frorth-DR6) = (1.464,0.315) with scat- 
ter 0.058 dex. Using the IDLspec2d reductions in- 
stead gives (adirect-IDLspec2d, ^direct-IDLspec2d) = 

(1.122,0.304), with intrinsic scatter 0.100 dex, and 

(a rth-IDLspec2d, ^orth-IDLspec2d) = (1.401,0.317) with 

scatter 0.061 dex, respectively. This suggests that we should 



assume as typical systematic error 5a ~ 0.05 and 5b ~ 0.02. 
These values are larger than the random errors quoted 
above, showing that is important to separate random from 
systematic errors. 

2.2 Dependence on selection 

If we do not account for selection effects we find mean val- 
ues (Mr) = -21.94, (log R e ) = 0.62, (fi e ) = 19.71 and 
(log a) — 2.30, and covariances Cu = 0.2660, Crr = 
0.0488, Cvv = 0.0127, Cm = 0.0820, Civ = 0.0036 and 
Crv — 0.0159. Although these covariances have changed, 
the coefficients of the Fundamental Plane change little: 
(^direct, ^direct) = (1.1723,0.2924) with intrinsic scatter 
0.079, and (a or th, 6 or th) = (1.4235,0.2914) with scatter 
0.047. Note that the intrinsic scatter is smaller! This is 
largely because the distribution of luminosities in a mag- 
nitude limited survey is narrower than in a volume limited 
survey (because the faint objects are under-represented); if 
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Figure 4. Fundamental Plane residuals with respect to the orthogonal fit in the g, r, i and z bands, shown as a function of distance 
Xpp along the Plane. Symbols show the median residual in narrow bins in Xpp; dashed and dotted lines enclose 68% and 95% of the 
points. Solid lines show the expected trend due to correlated measurement errors. 



one does not account for this, then the magnitude limited 
catalog is more homogeneous, making for a tighter FP. The 
fact that the coefficients a and b hardly change is largely 
fortuitous (Sheth & Bernardi 2009 show why they are ex- 
pected to change); it happens because the scaling relations 
are actually slightly curved (see Section 2.4, and also Hyde 
k Bernardi 2009). 

Had we not removed objects with b/a < 0.6 then 
Cu = 0.3756, C RR = 0.0549, Cvv = 0.0171, Cm = 0.0948, 
Civ = -0.0132 and C RV = 0.0171. The resulting FP 
has (adirect, ^direct) = (1.1674,0.2936) with intrinsic scat- 
ter 0.091, and (a or th, b ort h) = (1.4294,0.3052) with scatter 
0.055. Since objects with b/a < 0.6 are almost certainly not 
early- types, the fact that a and b are almost the same as 
when these objects have been removed suggests that the 
presence of a few non-early-types in our sample has little 
effect on our findings. 

Finally, selecting only galaxies with SDSS-DR6 velocity 
dispersions (i.e. excluding objects with low S/N spectra or 
presence of weak emission lines, i.e. the status flag not- 



equal to 4) changes a and b very slightly - the variation is 
less than 0.01 for both coefficients. Table 1 compares the r- 
band FP coefficients associated with these various selection 
and parameter choices. 

2.3 Dependence on waveband 

The Fundamental Planes associated with the orthogonal fits 
to the g, r, i and z band data are shown in Figure 3. The 
coefficients (a, b) of these Fundamental Planes are shown in 
each panel; they are also reported in Table 2, as are the cor- 
responding coefficients of the direct fits. Notice that there 
is a small but systematic increase of a with wavelength. To 
estimate its significance, we require an estimate of the errors 
on a and b. Although systematics associated with how a was 
measured make 5a ^ 0.05, 5b ^ 0.02, this additional system- 
atic error is not relevant if we wish to compare the Planes 
at different wavelengths, because the same choice for a is 
made for all wavelengths. What matters here is the typical 
uncertainty due to random errors on these best fit values. 



Stellar Mass Fundamental Plane 7 



Table 1. Dependence of r-band Fundamental Plane coefficients 
on sample selection and parameters. Typical uncertainties on the 
coefficients due to random errors are Sa ~ 0.02 and Sb ~ 0.01; 
sytematics are Sa ~ 0.05 and Sb ~ 0.02. 



Band 


a 


b 


rms mt 


direct 








r 


1.1701 


0.3029 


0.0964 


DR6-cr 


1.1892 


0.3032 


0.1081 


spec2d-cr 


1.1223 


0.3041 


0.1002 


no-Vmax 


1.1703 


0.3036 


0.0792 


all b/a 


1.1674 


0.2936 


0.0913 


orthog 








r 


1.4335 


0.3150 


0.0578 


DR6cr 


1.4642 


0.3151 


0.0581 


spec2dcr 


1.4013 


0.3173 


0.0613 


no-Vmax 


1.4235 


0.2914 


0.0473 


all b/a 


1.4294 


0.3052 


0.0554 



Table 2. Coefficients (a, b) of the Luminosity Fundamental Plane. 
Typical uncertainties on the coefficients due to random errors are 
Sa ~ 0.02 and Sb ~ 0.01. Sytematics errors give Sa ~ 0.05 and 
Sb ~ 0.02. 



Band 



rms obs rmsint 



direct 



9 


1.1154 


0.2957 


-8.0463 


0.1102 


0.1005 


r 


1.1701 


0.3029 


-8.0858 


0.1074 


0.0964 


i 


1.1990 


0.3036 


-8.0481 


0.1067 


0.0950 


z 


1.2340 


0.3139 


-8.2161 


0.1052 


0.0921 


orthog 
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1.4043 


0.3045 


-8.8579 


0.0696 


0.0617 


r 


1.4335 


0.3150 


-8.8979 


0.0664 


0.0578 


i 


1.4572 


0.3182 


-8.8914 


0.0652 


0.0563 


z 


1.4735 


0.3295 


-9.0323 


0.0635 


0.0538 



Section 2.1.2 and Figure 1 shows that these random errors 
are 5a ~ 0.02, 5b ~ 0.01. Thus, we have a 3a detection of 
the steepening of the Fundamental Plane with wavelength. 



2.4 Evidence for curvature 



Figure 4 shows residuals 

log R e 



a log a — b/j e — c 



(13) 



VI + o? + b 2 
from the orthogonal fit as a function of distance 

_ alogRe + logcr + biie/a 

A F p = (14) 

Vl + a 2 

along the Plane. Weak trends are seen in all bands. 

Interpretation of these trends is complicated by the fact 
that the measurement errors are correlated, and they are 
larger at small sizes. To see what effect is expected, we must 
use the numbers from the covariance matrices T and 6 to 
estimate (AfpXfp)/(X| p ), the expected slope of the cor- 
relation shown in Figure 4. In the r-band, (Afp Afp) = 
— 0.0007, where the first term is the contribution from T 
and the second from S. I.e., there is no intrinsic correla- 
tion (essentially, by definition), but there is a contribution 



from the measurement errors. This contribution is easy to 
estimate because Afp is almost proportional to IP — aV, 
where IP = \og 10 (R e ) — 0.3/x e , and the error in IP is well- 
known to be negligible, at least for deVaucouleur-like profiles 
(Saglia et al. 1997; or see bottom panel of Figure 5 in Hyde 
& Bernardi 2009), so the error in Afp is dominated by the 
error in V. On the other hand, the error in Afp is due to 
errors in R, V and IP. Errors in R and V are almost uncor- 
rected, so the dominant contribution from correlated errors 
comes from the fact that V appears with different signs in 
Afp and Afp- 

A similar analysis of (ApV) shows that it is dominated 
by the contribution from T rather than S. In the r-band, 
(Xfp) = 0.087 + 0.008. Thus, correlated errors are expected 
to produce a weak trend, with slope —0.007. The solid line in 
the Figure shows this trend - it is similar to that observed, 
suggesting that correlated errors can account for most of 
the observed weak decline in Afp. On the other hand, cor- 
related errors cannot account for the break downwards at 
small Xfp: this is genuine curvature. 

A similar analysis in the other bands yields the solid 
lines shown in Figure 4. This illustrates that sample sizes are 
now large enough that correlated errors produce systematic 
effects which must be accounted for when performing the 
fits. 

At the low size/mass end, the Plane is probably warped 
(Figures 4, and see also Figure 10 below). However, we view 
this with caution since uncertainties/systematics in the ve- 
locity dispersion/size measurements are also larger. The pos- 
sibility of curvature induced by the presence of spirals in the 
sample is discussed in Appendix A. We argue that the loca- 
tion of spirals on the fundamental Plane would not result in 
negative residuals at small Xfp , and that the spirals are to 
few to significantly bias our results. 



2.5 Variation in thickness along the Plane 

Notice that the Plane is thinner at large Xfp . Figure 5 shows 
that the typical measurement errors vary little with dis- 
tance along the Plane. Since the observed scatter is larger at 
smaller Afp (the width of the regions enclosed by the dashed 
and dotted curves in Figure 4 increases at small Afp), but 
the contribution from errors is approximately constant, we 
conclude that the intrinsic scatter around the Plane de- 
creases dramatically as Afp increases. It is remarkable that, 
at Afp > 1 , measurement errors account for essentially all of 
the observed scatter (about 0.04 dex), suggesting that the 
Plane is rather thin at large Afp. 



2.6 Potential biases from cuts in L 

There has been recent interest in the fact that the coeffi- 
cients of the FP depend on how the sample was selected, a 
point emphasized by Bernardi et al. (2003). E.g., D'Onofrio 
et al. (2008) and Nigoche-Netro et al. (2008) show that a 
decreases if faint galaxies are removed from the sample, and 
one does not account for the fact that they are missing. 
Donofrio et al. suggest that this may reflect the fact that lu- 
minous and faint galaxies had different formation histories. 
We show below that, although the latter may be true, the 
dependence of a on sample selection alone is not a reliable 
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Figure 5. Contribution of measurement errors to the residuals 
from the orthogonal fit as a function of Xpp for the r-band (sym- 
bols and line styles same as previous figure). Results for other 
bands are similar. 



indicator. We do so by constructing an FP relation using 
pairwise correlations with no curvature - so one would have 
concluded that nothing special was happening to galaxies 
at either end of the sequence. We then show that cuts in L 
change a and b in ways that are quite similar to that ob- 
served by D'Onofrio et al. In the following subsection, we 
show that similar effects also occur if cuts in a are made. 

Before discussing the FP, it is easier to first consider the 
size-surface brightness relation. At fixed L, (log 10 R\ji e ) will 
be a line of slope 1/5 with no scatter, by definition. Changing 
L moves this line to the left or right, but does not change 
the slope. However, because there is a correlation between 
R e and L, more luminous galaxies are bigger on average, the 
high L galaxies only populate the large R part of their line; 
galaxies of lower L only populate the lower part of their line. 
Thus, in general, the full log R — \i e correlation, which is got 
by averaging over the full range in L, will have a different 
slope than 1/5, with the difference being determined by the 
strength of the log R e — log L correlation. 

Figure 6 shows this explicitly using three narrow non- 
overlapping bins in L. In any one bin, the slope of the rela- 
tion (shown by the solid lines) is 0.2. The slope obtained 
from combining two nearby bins will be slightly steeper; 
combining all three bins would yield an even steeper slope. 
This is a generic argument, but note that it also works when 
(\ogR e \ logL) increases linearly with logL. Since the slope 
of the log R e — /i e correlation is 1/5 for a narrow bin in L, 
but something else when all L are included, the slope of the 
log R e — ii correlation will appear to depend on the range 
of L included in the sample even though the fundamental 
underlying correlation is linear. This, essentially, is the ori- 
gin of the behaviour seen by D'Onofrio et al. (2008) and 
Nigoche-Netro et al. (2008). 

To show this explicitly, the upper of the two dot-dashed 
lines in the Figure shows the result of averaging over the 
objects in the three luminosity bins shown in Figure 6: it is 
steeper than the relation for any one of the bins. Of course, 
in a magnitude limited survey, one does not simply aver- 
age over all the luminosity bins. Rather, when objects with 
the full range of L are used, then each object is weighted 



o 




19 20 ^ 21 

Me [mag/arcsec 2 ] 



22 



Figure 6. Effective radius versus surface-brightness in magnitude 
bins of (from botton to top) -20.5 < M r < -20, -22 < M r < 
-21.5 and -23.5 < M r < -23. The solid lines are linear fit to the 
galaxies in the different bins. The slope is 0.2 in all cases. Dashed 
and dotted lines show the regions which contain 68% and 95% of 
the objects in each bin. Upper dot-dashed curve shows a fit to the 
objects in all three bins, and lower dot-dashed curve shows the 
correlation in the full data-set which accounts for selection effect 
(as in Table 1 of Hyde & Bernardi 2009). 



by V~^(L). The lower dot-dashed line shows this relation: 
this shows that both the slope and the zero-point of this 
correlation are sensitive to this weighting. 

Extending this argument to the Fundamental Plane is 
slightly more involved. Once again, at fixed L, the scaling 
between R and \i is straightforward. The inclusion of a a- 
dependent term to the x-axis serves to shift the lines as- 
sociated with different a horizontally - the goal of the FP 
algorithm is to shift them so they lie on top of one another 
as much as possible. To see the effect this has, note that for 
a given L, each object with size R e is shifted by an amount 
which depends on its a, so the full Fundamental Plane con- 
sists of shifting each point horizontally in the log R e — \i e 
Plane by (log a\ log L, log R e ) on average. If the mean shift 
depends on R e , this will change the slope of the line of fixed 
L from 1/5 to something else. If there is scatter around this 
mean shift, then the line of fixed L will be broadened. The 
combination of change in slope and additional scatter both 
serve to make the parameters a and b of the Fundamental 
Plane depend on the range of L which are included in the fit, 
even though the underlying pairwise correlations are linear 
(i.e., their slope does not depend on L). 

Figure 7 shows all this explicitly. The symbols show 
how the FP coefficients change as faint objects are excluded 
(symbols on the left) or as luminous objects are excluded 
(symbols on the right). Large symbols show the effect of 
accounting for selection effects (each object is weighted by 
the inverse of the volume over which it could have been 
observed, so luminous objects have smaller weights), and 
smaller symbols show results when all objects are weighted 
equally (selection effects are ignored). Note that a depends 
strongly on the luminosity cut, whereas b is less strongly 
affected. The smooth curves show the result of making sim- 
ilar measurements in a mock catalog constructed following 
methods given in Bernardi et al. (2003). Note in particular 
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Figure 7. Dependence of FP coefficients on the luminosity range 
of the sample. The set of symbols on the left show how a and 
b change if low luminosity objects are excluded from the sample 
(the symbols show the absolute magnitude of the faintest object 
which is kept in the sample); the symbols on the right show what 
happens if high luminosity objects are excluded (symbols show 
the absolute magnitude of the brightest object in the sample). 
Large symbols show the result of accounting for the selection ef- 
fect which comes from the apparent magnitude limit of the SDSS; 
smaller symbols show the result of ignoring this effect. Smooth 
curves show similar measurements in a mock catalog in which all 
underlying correlations were pure power-laws. 

that all scaling relations in the mock catalog were linear; 
there was no curvature. Nevertheless, Figure 7 shows that 
the FP coefficients in the mock catalog depend on the value 
of the luminosity threshold similarly to how they do in the 
data. This demonstrates that a detection of dependence of a 
on luminosity threshold does not imply that the underlying 
scaling relations are curved. 

The strong dependence of a on luminosity threshold is 
easily understood. At fixed L, R and a are ant i- correlated; 
it is only when averaged over a large range in L that R and 
a are positively correlated (e.g. Bernardi et al. 2003). Since 
a is essentially the slope of the R — a relation, the result 
of restricting the range in L is to drive a to smaller values, 
since at fixed L it must be negative. This also explains why, 
when no account is taken of the magnitude limited selection, 
a changes little (compared to when the selection effect is 
accounted for) when low luminosity galaxies are excluded, 
but more strongly when high luminosities are excluded. In 
a magnitude limited sample, the luminosity distribution is 
biased to larger L (since they can be seen to larger volumes) , 
so the effect of removing low L is less severe (this removes a 
smaller fraction of galaxies), but the effect of a cut at high 
L is more dramatic (this removes a larger fraction of the 
galaxies) . 



Figure 8. Dependence of FP coefficients on the range of velocity 
dispersions in the sample. The set of symbols on the left show 
measurements in samples in which objects with a less than the 
value indicated were excluded; the set of symbols on the right 
show results when the sample excludes objects with a larger than 
the indicated value. Large symbols show the result of accounting 
for the selection effect which comes from the apparent magnitude 
limit of the SDSS; smaller symbols show the result of ignoring 
this effect. Smooth curves show similar measurements in a mock 
catalog in which all underlying correlations were pure power-laws. 

samples are selected on the basis of absolute magnitude, but 
many samples exclude objects with small a, simply because 
small velocity dispersions are difficult to measure. 

Figure 8 shows that b is (again) hardly affected. How- 
ever, a or th depends strongly on this selection, and the trend 
is opposite to that for luminosity, while afreet is less affected 
(also see Bernardi et al. 2003; Nigoche-Netro et al. 2008). 
This is easy to understand (also see discussion in Bernardi 
et al. 2003): afreet is close to the slope Crv/Cvv of the 
(logi?| log a) relation. If this relation is linear, then exclud- 
ing small or large values of a should not change this slope. 
On the other hand, a or th measures the orthogonal slope of 
this relation, and this must steepen as the bin in a narrows; 
after all, in the limit of a very narrow bin in cr, this orthog- 
onal slope must become infinite. (Sheth & Bernardi 2009 
show this analytically - a or th depends both on Crv/Cvv 
and on Cvv itself. Reducing the range of a reduces Cvv, 
thus increasing a or th-) Once again, similar measurements in 
the pure power- law mock produce similar trends, although 
in this case the agreement with the SDSS is not as good. 
Most of this (small) discrepancy can be attributed to the 
fact that the underlying pairwise scaling relations in the 
SDSS are slightly curved (see next section), an effect which 
we have deliberately removed from our mocks so as to il- 
lustrate that Figures 7 and 8 are not good diagnostics of 
curvature. 



2.7 Potential biases from cuts in a 

There is no fundamental reason why we should have re- 
stricted our study to how a and b depend on L. This sub- 
section shows how the FP coefficients change as the range 
of a in the sample is varied. This study is potentially more 
interesting than that of the previous subsection, since few 



2.8 Comparison with previous work 

Our fits differ slightly from those reported by Bernardi et al. 
(2003). This is not unexpected, because the SDSS improved 
its photometric reductions significantly between DR1 (on 
which Bernardi et al. 2003 was based) and DR6 (on which 
the present analysis is based). In addition, Figure 8 shows 
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Figure 9. Comparison of the stellar mass estimates from spectra 
(Gallazzi et al. 2005) with those from broad-band photometry 
(Bell et al. 2003). Symbols and line-styles same as in Figure 4. 

that simply removing a < 100 km s _1 increases a or th by 
about 15%. This is interesting because Bernardi et al. (2003) 
removed objects with a < 90 km s _1 from their sample (due 
to the dispersion of the SDSS spectrograph). The distribu- 
tion of velocity dispersions is relatively narrow, so this cut 
can have a non- negligible effect. Figure 8 shows that this 
might bias a high by about 10%, which is about the level 
of discrepancy between the Bernardi et al. (2003) FP and 
this work (but note that the maximum likelihood method of 
Bernardi et al. did try to account for the cut in a.) 

While the absolute values of the best-fit coefficients have 
changed, the relative values have not changed significantly: 
both a and b increase while the intrinsic scatter decreases in 
the redder bands. This is true whether the fit minimizes the 
scatter in the log R e direction, or in the direction orthogonal 
to the Plane. In addition, it is thinner at the large \ogR e 
end. The weak but significant increase of a with wavelength 
is consistent with that reported by Bernardi et al. (2003) 
and, more recently, by La Barbera et al. (2008). 



3 FP*: THE STELLAR MASS FUNDAMENTAL 
PLANE 

The previous section studied the traditional Fundamental 
Plane, FP. In this section, we show the result of replacing 
luminosities with stellar masses to produce what we will call 
FP*. To illustrate that our results are robust to changes in 
how one estimates the stellar mass, we show results based 
on two different estimates of M*/L: one, from Gallazzi et 
al. (2005), which is based on a likelihood analysis of the 
spectra; another, from Bell et al. (2003), who use the k- 
corrected broad-band color (M*/L r ) = 1.097 (g-r) -0.306. 
Figure 9 compares these estimates with each other; there is 
an offset of about 0.15 dex and a small trend which shows 
that Moaiiazzi/MBeii tends to increase at larger stellar mass. 
Most of this offset (0.11 dex) is due to the difference in initial 
mass function (IMF) used in the stellar population models 



Table 3. Coefficients (a, (3) of the Stellar Mass Fundamental 
Plane. Typical uncertainties on the coefficients due to random er- 
rors are 5a ~ 0.02 and 5b ~ 0.01. Sytematics errors give 5a ~ 0.05 
and 5b ~ 0.02. 



M*/L 


a 


P 


c 


rms obs 




direct 












Spectra 


1.3989 


0.3164 


4.4858 


0.1160 


0.0894 


Color 


1.3501 


0.3293 


5.0015 


0.1115 


0.0835 


orthog 
Spectra 
Color 


1.6287 
1.5462 


0.3359 
0.3449 


4.4238 
4.9300 


0.0648 
0.0638 


0.0486 
0.0466 



by the two groups. (Gallazzi et al. used Bruzual & Chariot 
2003 models which assume the Chabrier 2003 IMF, while 
Bell et al. assume a diet-Saltpeter IMF.) The rms scatter 
between these two estimates increases significantly at small 
M*. While the offset will affect the zero-point but not the 
slope of the FP*, the weak trend could actually introduce a 
small systematic effect in the slope. 

Figure 10 shows the associated Fundamental Planes, 
the coefficients of which are reported in Table 3. The uncer- 
tainties in 5a and 5b are similar to those described in Sec- 
tion 2.1.2 (i.e., random errors 5a ~ 0.02 and 5b ~ 0.01, with 
larger systematic uncertainties of 5a ~ 0.05 and 5b ~ 0.02). 
In this case also, understanding the errors is important. We 
have assumed that the error in /x* is given by the error in 
the photometric quantity \i plus a contribution from the er- 
ror in logM*/L. Gallazzi et al. (2005) report rms errors of 
about 0.06 dex, so we add 2.5 times 0.06 dex in quadrature. 
We further assume that the error in logM*/L is uncorre- 
cted with that in the size or velocity dispersion estimates. 
Whereas the first assumption is probably correct, the second 
is almost certainly incorrect for the stellar mass estimates 
which come from the spectra. (Gallazzi et al. 2005 do not 
show how the error on M*/L correlates with the error on 
other observables. They do show that the M*/L error does 
not correlate with spectral type.) This may explain some of 
the differences we see between the two stellar mass Planes. 

The important point is that in both cases, the slope of 
FP* is steeper than that of the FP in the z-band, but it 
is shallower than the virial scaling. This suggests that the 
ratio of dynamical to stellar mass, M^ yn /M*, varies system- 
atically across the population. Our results suggest that it 
is a weakly increasing function of Md yn or M*. Figure 11 
shows a direct comparison: {M dyn /M*\M dyn ) oc M^ 7±0 - 01 
where Mdyn = 5R e cr 2 /G and M* is from Gallazzi et al. (the 
error 0.01 on the slope was computed accounting for system- 
atics errors - the uncertainty from random errors is smaller 
~ 0.003). If there were no scatter around this relation, then 
we would expect {M dyn /M* |M*) oc M* 017/0 - 83 oc Mf; be- 
cause there is scatter, this scaling is shallower, oc m®' 06±0 ' 01 
(Hyde & Bernardi 2009). Gallazzi et al. (2006) also measured 
the Mdyn, M* relationship using deVaucouleur radii of SDSS 
galaxies. They obtain (M dyn /M*\M dyn ) oc M^ ±om which 
is statistically equivalent to our result. 



Stellar Mass Fundamental Plane 11 




Figure 10. The stellar mass Fundamental Plane associated with the two esimates of the stellar mass (from spectra, left and photometry, 
right). Symbols and line-styles same as for Figures 3 and 4. Typical uncertainties on the coefficients due to random errors are 5a ~ 0.02 
and 5b ~ 0.01. Sytematics errors give 5a ~ 0.05 and 5b ~ 0.02. 



4 THE PLANE IN ^-SPACE 

Bender, Burstein Sz Faber (1992) suggested that the combi- 
nation of variables which makes up the Fundamental Plane 
could be combined in a more physically transparent way. 
They called this combination k- space. The three axes of this 
space are 



Kl = 



K3 = 



\og 10 (R e a 2 ) _ log 10 (a 2 / e 2 /i?e) 



V2 ' 
log 10 (i? e q 2 /L) 



K2 



V6 



, and 



(15) 



It is interesting to examine how early types are distributed 
in what we call k*- space, which is obtained by replacing 
luminosity with (Gallazzi et al. 2005) stellar mass in the 
expressions above. 

Figure 12 compares the two most interesting projections 
of k- (left) and ft*-space (right), when all variables have been 
expressed in solar units. The top panels show the 'edge-on' 



view: M dyn /L- M dyn and M dyn /M* - M dyn . Note that the 
top right panel is essentially the same as Figure 11, except 
that now we show the distribution using contours (objects 
have been weighted by 1/Vmax to account for selection ef- 
fects); dashed lines show forward and inverse fits, and the 
solid line, which has slope 0.383 and 0.429 in the two panels, 
shows the bisector fit. 

The bottom panels show the distribution of objects 
within the Plane. The dashed line in the panel on the left 
shows fti + K2 ^constant (as suggested by Bender et al. 
1992). In the solar units used in the plot, the dashed line 
shows 
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(16) 
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Figure 11. Ratio of dynamical and stellar masses (from Gallazzi 
et al. 2005) as a function of dynamical mass. Symbols and line- 
styles same as Figure 4. The thick solid curve is a quadratic fit 
(see Hyde & Bernardi 2009). 

whereas it is 

( M,/Rl \ (M dyn y /3 _ ( M dyn y 1/V ~ 3 
VlO lo M /kpc 2 / \ M r J ' \10 10 MqJ 

(17) 

for the panel on the right. In the case of stellar masses, it 
is helpful to cube both sides of this expression, and then 
rearrange so that all the dependence on the dynamical mass 
is on the same side. This yields 

U lo M /kpc 3 J U™oJ ' (18) 

evidently, the dashed line expresses a relation between stellar 
density and dynamical mass: the upper limit to the stellar 
density is approximately proportional to M7^ 3 . 



5 DISCUSSION AND CONCLUSIONS 

We showed that the slope of the traditional, luminosity- 
based Fundamental Plane FP depends weakly but systemat- 
ically on waveband - it steepens towards the virial relation 
in the redder bands (Figure 3 and Table 2). Replacing L 
with stellar mass M* leads to coefficients which are slightly 
closer to the virial ones; this is true whether one estimates 
M*/L from the spectra or simply from broad band colors 
(Table 3 and Figure 10). 

The stellar mass FP is slightly thinner (with an aver- 
age intrinsic scatter of ~ 0.048 dex) than the luminosity- 
based FP (which has a typical intrinsic scatter larger than 
~ 0.055 dex). The intrinsic scatter also decreases with wave- 
lengths: from 0.062 dex in g to 0.054 dex in z. 

We also showed that the intrinsic scatter around both 
the luminosity-based and stellar mass Planes becomes sig- 
nificantly broader at low sizes/masses (Figures 4 and 10) 
and that measurement errors account for essentially all of 
the observed scatter (about 0.04 dex) at large sizes/masses, 
suggesting that the Plane is rather thin for the very massive 
galaxies. 

The fundamental nature of the stellar mass Plane FP*, 



and the fact that g — r color is a good indicator of M*/L 
(e.g. Bell et al. 2003), explains why residuals from FP cor- 
relate with color (e.g. Bernardi et al. 2003), or, equivalently, 
why color may be thought of as the fourth Fundamental 
parameter in early- type galaxy scaling relations. 

The fact that FP* does not quite have the virial scal- 
ings suggests that the ratio of stellar to dynamical mass, 
M* / 'Mdyn should vary across the population. Figure 11 
showed this was indeed the case: Md yn /M* oc Mj^J. This 
is in qualitative agreement with the results of the SAURON 
project (Cappellari et al. 2007), which is based on a very 
different analysis technique. At higher redshifts, the Funda- 
mental Plane method is technically less challenging, so we 
expect it to provide a useful measure of the evolution of 
Mdyn I M*. This will also provide a useful check on the sug- 
gestion that the slope of the correlation between Md yn and 
M* does not evolve out to z ~ 1 (Bundy et al. 2007). 

We also presented an analysis of ft-space (Bender 
et al. 1992), but after replacing luminosities with stellar 
masses (using values from Gallazzi et al. 2005). The plot 
of M* I Mdyn versus Md yn referred to above is also known 
as the edge-on view of ft* -space. The face-on view of the 
Plane in this space (Figure 12) showed that the maxi- 
mum stellar density is smaller in the more massive objects: 
M*/i?f oc M^y^ 2, . The ft*-space scalings at low redshift, and 
an estimate of how they evolve, should provide interesting 
constraints on galaxy formation models. 

Datasets are now sufficiently large that statistically sig- 
nificant curvature in most scaling relations has now been 
seen (e.g. Hyde & Bernardi 2009). However, illustrating that 
the Fundamental Plane is warped is more difficult. Previous 
work has shown that the coefficients of the FP can change 
dramatically as objects of different type (e.g. low luminosity 
or a) are removed from the sample (e.g. Bernardi et al. 2003 
for dependence on a, and Donofrio et al. 2008 and Nigoche- 
Netro et al. 2008 for dependence on L). This has lead to 
the suggestion that these changes indicate that the FP is 
warped (e.g. Donofrio et al. 2008). We showed that similar 
changes arise even in samples where the underlying pairwise 
scaling relations are not curved - i.e., when the Plane is not 
warped (Sections 2.6 and 2.7). Thus, the dependence of the 
fit parameters on the range of L or a in the sample is not, by 
itself, evidence for curvature. This dependence on the range 
of L or a may explain some of the relatively wide range of 
Fundamental Plane coefficients in the literature. 

A more robust measure of how warped the Plane was 
also discussed. We showed that a good understanding of the 
errors is necessary to interpret the results; else, correlated er- 
rors might lead one to conclude there is curvature even when 
there is none (Figure 4 and discussion in Section 2.4). We 
conclude that the Plane is warped at the low size/mass end 
(Figures 4 and 10), and that it is also significantly broader 
at low sizes/masses. These conclusions are not affected by 
a low rate of spiral galaxy contamination, as discussed in 
Appendix A. 

Our analysis of the bias introduced in the FP coeffi- 
cients by cuts in L and a raises an interesting puzzle. The 
~ 50 objects in the SLACS sample (Bolton et al. 2008) were 
selected because they are gravitational lenses, meaning that 
they tend to have large velocity dispersions. The sample is 
particularly interesting because dynamic, stellar and gravi- 
tational lensing mass estimates are available. Although the 



Stellar Mass Fundamental Plane 13 




10.0 10.5 11.0 11.5 
Log 10 M dyn [M Sun ] 



12.0 



9.5 



9.0 



S. 8.5 



o 




1.2 
1.0 
0.8 
\ 0.6 

■o 

° 0.4 
cn 
o 
_l 

0.2 
0.0 
-0.2 



9.5 



9.0 



8.5 



I I I I I I I I I 

- 


i i i i | i i , i | i i i i | i iv i i 

/ 

/ 

/ 

/ ~ 




/ 

/ /- 
/ / . 

/ / ~ 

/ / 
/ / 




■ / / - 

! l \ ;' 4 ••. ••. •-. •• •". '• • / ; ; : ; ' . • ■ ; • 

. '. .' ' — Z 
.:-.-. ^*s^ ^ ^ - - 










pfT: • .•' ••' .y/- ' ' '' 

'»,.: ..• . y.. .• / 


/ /^••- '•.*.'• •.'.'•:•' - 


"::::::. ' ' '•' 1 i i !• ;J/T . • ■ ■ 
1 ! ' f ■ 




...XI../.. 


. . . . 1 . . , . 1 . . . . 1 . . . . 


10.0 10.5 11.0 11.5 12.0 


Log 10 M dyn [M Sun ] 




Ill 

s 

N 

N 

.:•".•*'•.. N 

'."•":••.'.*■'•'••'« n 
*. • ■ ' " • n. ■ 




."'••s - - 

• — ~ 

/ /" — ' — '•'•*•' •- 

^yh' X- 








k > — 1 — ' ' :.• y : n . 






i i i i I i i i i 


i i i i 1 i i , i 1 i i i i 1 i i i i 



10.0 10.5 11.0 11.5 12.0 
Log 10 M dyn [M Sun ] 



10.0 10.5 11.0 11.5 
Log 10 M dyn [M Sun ] 



12.0 



Figure 12. The tt-space Fundamental Plane in r band (left panels) and with stellar masses (from Gallazzi et al. 2005) in place of 
luminosity (right panels). The figure shows «i, K2 and ^3 transformed to solar units. Top panels show versus dashed lines 

show the direct and inverse linear fits and solid line the bisector one. Bottom panels show versus v^^i- Dashed line sloping 

down and to right shows n\ + ^2 = constant scaled to solar units. 



small sample size means the intrinsic thickness of the FP in 
this sample is almost certainly underestimated (Figure 2), 
what is of interest here is the fact that the smallest velocity 
dispersion in this sample is a ~ 160 km s _1 . 

Bolton et al. make no attempt to correct for the fact 
that their sample has no objects with small a. Instead, they 
argue that, given the same size and velocity dispersion distri- 
butions, the SLACS sample has approximately the same lu- 
minosity distribution as the parent SDSS sample from which 
it was drawn. They also report that the coefficients of the 
Fundamental Plane for their sample are rather close to the 
virial scalings, and they use this to motivate a number of 
conclusions about the origin of 'tilt' in most early- type sam- 
ples. However, our Figure 8 suggests that if objects with 
a < 100km s _1 are missing from the sample (this is a fairly 
standard cut-off based on what most spectrographs used for 
this work are capable of), and no account is made of this, 



then a or th is biased high by about 15% (also see discussion 
in Section 2.2 as well as Figure 4 in Bernardi et al. 2003). If 
objects with a < 160 km s _1 are missing, then a or th is bi- 
ased to values which are close to 2! Thus it is possible that 
the slope a or th in the SLACS sample is high simply because 
it is biased. If so, then the conclusions about the origin of 
'tilt' should be revised. 
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APPENDIX A: SAMPLE SELECTION AND 
DISK CONTAMINATION 

The main text describes how our sample was selected; the 
main cut is on the shape of the light profile fracDev = 1. 
However, a non-negligible fraction of these objects (about 
20%) have b/a < 0.6, which we remove because we believe 
the axis ratio is caused by an edge-on disk component. Pre- 
sumably, similar objects viewed faced-on remain in our sam- 
ple, so the question arises as to how they may have affected 
our results. 

To address this, we divided the full sample in 5 bins in 
luminosity, randomly selected ~ 200 galaxies in each bin, 
and classified them, by eye, as E/S0, Spirals, and Others 
(irregulars, low surface brightness objects, interacting sys- 
tems, projections). The different panels in Figure Al show 
these classifications in four luminosity bins. In each panel, 
the solid histogram (the one with the most counts) con- 
tains rsj 200 objects (a few objects were removed because 
of contamination of a bright star or misclassification) . The 
dot-dashed histogram shows the subset with f racDev = 1 in 
both g and r bands, the dashed histogram shows the result 
of applying all our other cuts except the one on b/a, and the 
dotted shows the result of excluding objects with b/a < 0.6. 
This shows that f racDev = 1 removes most of the Spirals 
and Others - our further cuts help in reducing the left over 
contamination. 

While our selection removes most of the later type 
galaxies, it is perhaps surprising that a significant fraction 
of the objects we classified as ellipticals by eye do not have 
f racDev = 1. A selection of these objects is shown in Fig- 
ure A2; it is evident that many of these objects are only 
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Figure Al. Effect of selection cuts on the morphological composition in our sample, shown for four bins in absolute magnitude. Solid 
histogram shows everything, dot-dashed is the subset with fracDev = 1, dotted is the smaller subset which has spectra and satisfies our 
cuts in eClass, and dashed is for b/a > 0.6. 



marginally E/SOs, so we are confident that fracDev = 1 is a 
reliable cut. To illustrate the effect that such objects might 
have had on our results, the diamond-shaped symbols in Fig- 
ure A3 show their location relative to the FP denned in the 
main text. These objects are offset to higher sizes, and they 
have substantially higher scatter around the FP (see inset 
at top left), suggesting that our decision to exclude them 
from the analysis in the main text, is reasonable. 

The filled circles and open squares in Figure A3 show 
objects we classified as E/SOs and Spirals, and which satis- 
fied our selection cuts on fracDev and b/a (i.e., they make 
up the dotted histogram in Figure Al). There are a handfull 
of spirals and they tend to have larger R e - as one might 
expect from trying to fit a disk component with a single de- 
Vaucouleur profile. They are a slightly larger fraction of the 
total counts at small cr+0.21/i e , however they are so few that 



they do not bias the slope, scatter, or curvature properties of 
the FP. Additionally, the curvature described in Section 2.4 
results in negative residuals at low Xfp- Any curvature in- 
troduced by spiral contamination would introduce positive 
residuals because of artificially large radii. 



16 J. B. Hyde & M. Bernardi 



J1G221 2,76+432936.3 


JT2495fl.4a*fi42 103.2 


.1102746,41+112202,7 


J15132&34+025634,7 


J120-[23,92+t3183B 






E 






• 


• 


P 


J12D107,94+f32332.8 


JQ64746.96+522411.9 


J 125853, 78-0 308 137 


JQ95539,81 +39201 3,6 


J1fS27f 2,74+46361 8.3 




• 


♦ 




• 


J123801,99+4246G&2 


J 155636 86+415413,6 


J1 30^,33+650030.6 


J0B5S46, 12+041 11 1.5 


Jl 31 616.92-563457 .6 


• 


• 






• 


• 

• • 




J16O3m97+550645.3 


J1G2337, 12+1 25433,3 


J075141,64+293214 


J0B451 9*295327.4 


JOf7423JO,07+t94702.2 




• 


- 




• 


J16O01 7,06+273566,6 


J1Q3917,3+4745G8,9 


J 1 4540 1 .32- 6323 1 4.5 


J212735 83-065311 1 


J171717,fl7+301613.3 


mm 


• 




% 


• 

• 



Figure A2. Selection of objects classified by eye as E/SOs, but which have f racDev < 1. Objects such as these are shown by the diamond 
symbols in Figure A3. 
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Figure A3. The Fundamental Plane in the r band (same as top right panel in Figure 3). Filled circles and open squares show objects 
which satisfied our selection cuts and were morphologically classified by eye as E/SO (filled circles, see also dashed line in the inset at top 
left) and Spirals (open squares). Diamond-shaped symbols (and dotted line in the inset at top left) show objects we classified as E/SO 
by eye but do not have fracDev = 1. 
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